New analysis method of the halo phenomenon in finite many-fermion systems 
First applications to medium-mass atomic nuclei 
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A new analysis method to investigate halos in finite many-fermion systems is designed, as existing 
characterization methods are proven to be incomplete/inaccurate. A decomposition of the internal 
wave-function of the A^-body system in terms of overlap functions allows a model-independent 
analysis of medium-range and asymptotic properties of the internal one-body density. The existence 
of a spatially decorrelated region in the density profile is related to the existence of three typical 
energy scales in the excitation spectrum of the (A— l)-body system. A series of model-independent 
measures, taking the internal density as the only input, are introduced. The new measures allow a 
quantification of the potential halo in terms of the average number of fermions participating to it and 
of its impact on the system extension. Those new "halo factors" are validated through simulations 
and applied to results obtained through energy density functional calculations of medium-mass 
nuclei. Performing spherical Hartree-Fock-Bogoliubov calculations with state-of-the-art Skyrme 
plus pairing functionals, a collective halo is predicted in drip-line Cr isotopes, whereas no such effect 
is seen in Sn isotopes. 

PACS numbers: 21.10.Gv, 21.10.Pc, 21.60.Jz 
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I. INTRODUCTION 

The study of light nuclei at the limit of stability has 
been possible in the last two decades thanks to the 
first generations of radioactive ion beam facilities. One 
of the interesting phenomena observed close to the nu- 
cleon drip-line is the formation of nuclear halos. In such 
systems, either the proton or the neutron density dis- 
plays an unusually extended tail due to the presence of 
weakly-bound nucleons Since the first experimen- 
tal observation of such an exotic structure in ^^Li 0, 
other light neutron halo systems have been identified, e.g . 
6He i, "Be [5|,y,Q, "Be ii, ^^B i or i^C i^. 
On the proton-rich side, theoretical works demonstrated 
the existence of halo structures in spite of the presence 
of the Coulomb barrier 11 ill , as was seen experimen- 
tally for 8b [H [H [H, [ijand ^^Ne [H [Tj- Halos 
in excited states have been observed for ^^F [l^, 18 1, 
^^B [13 or "B [13, and several others are predicted O. 
It is worth noticing that weakly-bound systems ex- 
tending well beyond the classically-allowed region have 
also been theoretically predicted or experimentally ob- 
served for molecules (^He-^He-^^K [1|, ^Hea [24, 25, 26], 
•^IIe'*He2 [2j...). atom-positron complexes (e+Be, PsLi+, 
PsHe+...) [ill and hypernuclei (^H) f29l|. 
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The theoretical description of light halo systems is 
rather well under control. It usually relies on a cluster vi- 
sion where one (^^Be, ^^C.) or two (^^Li, ^He...) loosely 
bound nucleons define a low-density region surrounding 
a core. Assuming that core and halo degrees of freedom 
can be decoupled, essentially exact solutions of the sim- 
plified many-body problem are obtained by solving the 
Schrodinger equation for two-body systems Isol. [3l|, o r 
Faddeev equations for three-body ones [1, l30l l32l 
However, the boundary between halo and non-halo nu- 
clei is blurred by the presence of core excitations. Indeed, 
the inert decoupling of the loosely bound nucleons from 
the core is only an approximation. Nevertheless it has 
been assessed that halo systems arise when [34, 35] (i) 
the probability of nucleons to be in the forbidden region 
outside the classical turning point is greater than 50%, 
and (ii) the cluster structure is dominant and accounts 
for at least 50% of the configuration. Such conditions 
have been thoroughly studied [H, 113 and found to be 
fulfilled when (a) the separation energy of the nucleus is 
very small, of the order of 2 MeY/A'^^'^ , (b) the loosely 
bound nucleons occupy low angular-momentum states 
(£ = or £ = 1) for two-body clusters, or low hyperangu- 
lar momentum states (K = or if = 1) for three-body 
ones in order to limit the effect of the centrifugal barrier 
that prevents nucleons from spreading out [S^, and (c) 
the charge of the core remains small for proton halos. 
The latter requirement might be weakened because of a 
potential Coulomb-induced rearrangement of the single- 
particle states [39|]. 



2 



Going to heavier nuclei, few-body techniques face the- 
oretical and computational limits because of the large 
number of degrees of freedom involved. Single-reference 
energy density functional (SR-EDF) methods under the 
form of self-consistent Hartree-Fock-Bogoliubov (HFB) 
calculations become appropriate [ioL [4ll |. The EDF, ei- 
ther non-relativistic (Skyrme (42 . or Gogny 44j]) or 
relativistic [H, HI, [43, IH, E^, constitutes the only 
phenomenological input to the method. Phenomeno- 
logical functionals have now reached an accuracy suit- 
able for comparison of various observables with ex- 
perimental data over the known part of the nuclear 
chart [13, HH, m, HI]. However, properties of current 
EDFs are not yet under control in extreme conditions, 
where low-density configurations, isospin or surface ef- 
fects come strongly into play. Thus, the capacity of exist- 
ing functionals to predict properties of exotic nuclei, such 
as their limits of stability, remains rather weak f54!| . In 
that respect, the input from the coming generation of ra- 
dioactive beam facihties (FAIR at GSI, RIBF at RIKEN, 
REX-ISOLDE at CERN, SPIRAL2 at GANIL...) will 
help to further constrain models and to design a univer- 
sal EDF. 

Halo structures may contribute significantly to such 
a quest as they emphasize low-density configurations 
and surface/finite-size effects. Their study in medium- 
mass nuclei might provide relevant information regard- 
ing isovector density dependencies and gradient/finite- 
size corrections in the energy functional. In particular, 
the pairing strength in low density regimes and the evo- 
lution of shell structures towards the limit of stability 
might be constrained. However, two questions arise as we 
discuss potential medium-mass halos. Indeed, medium- 
mass nuclei are (i) large enough that the cluster picture at 
play in light nuclei needs to be revisited, in such a way 
that our understanding of the halo phenomenon might 
change significantly, and (ii) light enough that explicit 
correlations associated with symmetry restorations and 
other large amplitude motions are important and may 
impact halo properties. Including such correlations re- 
quire to perform multi-reference (MR) EDF calculations 
based on projection techniques and on the generator co- 
ordinate method (GCM) [iB, [H, [s^l . 

The first part of the present work, is dedicated to in- 
troducing a new method to identify and characterize halo 
structures in finite many-fermion systems. Although we 
only apply the method to even-even, spherical, medium- 
mass nuclei in the present paper, its range of applicability 
is wider [s^]. Regarding nuclei, extensions of the method 
to odd and deformed systems can be envisioned. The 
charge restriction for proton halos identified in light nu- 
clei is such that we do not expect proton halos in medium- 
mass systems. As a result, the present work focuses on 
exotic structures at the neutron drip-line. 

The article is organized as follows. Sec. |ll] provides 
a brief overview of the features that are crucial to the 
formation of halos. In Sec. IIIIl the limitations of exist- 
ing tools used to characterize skins and halos, such as the 



Helm model ^3], are highlighted. A new method to prop- 
erly identify and characterize halo features of weakly- 
bound systems in a model-independent fashion is intro- 
duced in Sec. IIVI We validate the method using a se- 
lection of toy models before applying it to the results of 
self-consistent spherical HFB calculations of Cr and Sn 
isotopes in Sec. |Vl The latter section is also devoted to 
a critical discussion of our results. Our conclusions are 
given in Sec. |Vll 



II. BASIC FEATURES OF HALO SYSTEMS 

The goal of the present section is to outline some of 
the elements that are crucial to the formation of halos. 
This will serve as an introduction to the more quanti- 
tative discussion proposed later on as we develop our 
new analysis method. For convenience, the discussion 
is conducted within the EDF framework whose basic as- 
pects are briefly recalled at first. Note however that the 
features discussed are not specific to a particular many- 
body method or approximation but constitute generic 
aspects of halos. For simplicity, spin and isospin indices 
are dropped in the present section. 



A. Elements of the nuclear EDF method 

The nuclear EDF approach is the microscopic tool of 
choice to study medium-mass and heavy nuclei in a sys- 
tematic manner [4l|. We consider a Single- Reference 
EDF formalism. In such an implementation, the energy 
is postulated under the form of a functional £[p,K,K*] 
of the (local or non-local) density p and pairing tensor 
K. The density matrix and the pairing tensor are further 
represented through a reference state 1$) 



($|c]c,|$) 
($1$) 



(1) 

(2) 



which takes the form of a quasiparticle vacuum and which 
reduces to a standard Slater determinant if no explicit k 
dependence of the EDF {£) is considered. Such a product 
state reads 



(3) 



where C is a complex normalization number whereas the 
quasiparticle operators {Pi, (3^) are obtained through the 
Bogoliubov transformation ([/, V) of the creation and an- 
nihilation operators {c\,Ci) defining an arbitrary single- 
particle basis 
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B. Importance of low angular- momentum orbits 



Pi = 



(4) 



The equations of motion, the so-called HFB equations, 
are obtained by minimizing the energy k, k*] with re- 



spect to the degrees of freedom (py , p. 



der the constraint that the neutron and proton numbers 
are fixed on the average in the reference state |<f>). This 
leads to solving the eigenvalue problem 



h-X A 

-A* -{h*-X) 



(5) 



where the one-body field h and the pairing field A are 
defined as 



d£ 



(6) 



A < being the chemical potential. Solutions of Eq. ([5|) 
are the quasiparticle eigenstates {U, V)^ whose occupa- 
tions are defined through the norm of the lower com- 
ponents Vv 



k 



\V,uV^ / |K(r)r dr 



(7) 



In order to analyze the properties of the many-body 
system, it is convenient to introduce the canonical 
hasis'^^^ {I'/'i)} HO) HOl- this basis, individual states 
can be grouped in conjugated pairs The one-body 

density p is diagonal whereas the pairing tensor k takes 
its canonical form 



Prj 



(8) 
(9) 



where Ui — > Q and Vi — —Vi play the role of BCS- 
like coefficients; vf being the canonical occupation num- 
ber. Even though the EDF method is not an independent 
particle theory, it is convenient to use the canonical ba- 
sis for analysis purposes as it provides the most intuitive 
single-particle picture and allows one to define individual 
"energies" and "pairing gaps" through 



We first discuss the impact of low-angular momentum 
orbitals*^^' on the density profile of halo nuclei. To do 
so, we first use the realization of the EDF method in 
which the reference state is taken as a Slater determi- 
nant. This corresponds to eliminating the dependence 
of the EDF on anomalous densities and thus the explicit 
treatment of pairing correlations. It is important to stress 
that, at least in principle, this does not mean that the 
effect of superfiuidity could not be accounted for in such 
a realization of the EDF method. It would, however, 
certainly require the design of more involve energy func- 
tionals that those used traditionally; i.e. Skyrme [i^ 
and Gogny [44] EDFs. 

Within such a realization of the EDF method, the 
HFB equations reduce to a standard one-body eigenvalue 
problem that provides the orbitals '^y{r) from which the 
auxiliary Slater determinant 1$) is constructed. Such a 
basis coincides in this case with both the canonical ba- 
sis and the quasiparticle basis introduced in Sec. Ill Al 
Restricting the description to spherical systems, consid- 
ering for simplicity a multiplicative local potential U{r) 
and forgetting about the spin degree of freedom, it can 
be proven [6l| that the density p{r) behaves asymp- 
totically as e"^'*" ""/(kq '')^) where the decay constant 
Kq = \/—2meQ/h? is related to the eigenenergy eo of the 
least bound occupied orbital in the reference Slater deter- 
minant. As the density used in the SR-EDF method is 
meant to reproduce the internal local density (see Ap- 
pendix lA 2\) . an analogue of Koopmans' theorem [62| 
holds, that is eo is equal to minus the one-nucleon separa- 
tion energy S'„ — Eq^^ ~ Eq , where Eq is the ground 
state internal energy of the A^-body system. As a result, 
long density tails arise for weakly-bound systems; i.e. in 
the hmit = |eo| — > 0. 

A more quantitative characterization of the density is 
provided by its radial moments (r"). Such moments are 
of special interest in the case of halo systems. At long 
distances, the dominant contribution to (r") comes from 
ifQ. In the limit of weak binding eg ~^ Oj the individ- 

ual moment (r")o (i) diverges as Cq ^ for n > 2£ — 1, 
(ii) diverges as In(ep) for n = 2£ — 1, or (iii) remains fi- 
nite for n < 2^ — 1 [61[. In particular, one finds that the 
wave function normalization (r'^)o diverges for s waves, 
whereas the second moment (r^)o diverges for both s and 
p waves. As a result, the root-mean-square (r.m.s.) ra- 
dius, defined as 



A,; 



(10) 

(11) 



^ The canonical basis is the name given to the natural basis in the 
context of HFB calculations. 



^ Although the notion of orbital often refers to an independent- 
particle picture or a Hartree-Fock approximation, it is important 
to note that the EDF method includes correlations beyond such 
approximations. In fact, and as discussed in Sec. IIVI the notion 
of orbital should rather be replaced by the one of overlap function 
in the present discussion. 
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diverges as eo — > if (po corresponds to a s or a p wave. It 

diverges as Eq ^ for a s wave and as Cq " for a p wave. The 
centrifugal barrier confines wave functions with higher 
orbital-angular momenta, in such a way that i?r.m.s. re- 
mains finite as eo ^ if ipo has an angular momentum 
£ >2. Equivalent arguments are found in the case of 
three-body systems (38j . 

According to the above analysis, only low- lying s or 
p waves near the threshold are able to extend signifi- 
cantly outside the classically forbidden region. The con- 
sequences of such patterns are that (i) one usually focuses 
on the evolution of the neutron r.m.s radius as a func- 
tion of neutron number, looking for a sudden increase 
as a signature of the building of a halo, (ii) the pres- 
ence and occupation of low-lying s or p waves are often 
seen as a prerequisite for the formation of neutron halos, 
(iii) orbitals with £ > 2 are not believed to contribute 
to halos. However, it is important to notice that (r^) is 
only the leading moment in the representation of the den- 
sity. The complete expansion of p{r) involves moments 
of higher orders which probe the nuclear density at in- 
creasing distances. Even if those higher-order moments 
weight usually little in the expansion, one cannot rule out 
£ >2-type halo structures, as (r")o with n > 2 diverges 
in the limit eg — > for such angular momenta: (r"*) di- 
verges for £ ^ 0,1, 2, (r^) diverges for ^ = 0, 1, 2, 3... and 
so on fssj . 



diverge, whatever the angular momentum of the least 
bound quasiparticle. In other words, the effect of pair- 
ing correlations is to induce a generic anti-halo effect by 
localizing the density [tO, [zH ■ 

Two additional effects may however blur such a pic- 
ture. First, recent HFB calculations performed in terms 
of a fixed one-body Wood-Saxon potential have shown 
that such a pairing anti-halo effect could be ineffective 
under extreme conditions [zl, [tI] . Indeed, very weakly 
bound Si/2 states (bound by a few keVs) tend to decouple 
from the pairing field because of their abnormal exten- 
sion. As a consequence, Eq = min^[i?i^] tends towards 
zero again as cq « A « and the r.m.s. radius of such 
an unpaired orbital may diverge, contributing strongly to 
the formation of a halo. Although this possibility is to be 
considered in principle, the depicted situation of a £ = 
orbit bound by a few keVs right at the drip-line is rather 
improbable and would be highly accidental in realistic nu- 
clei. Second, the pair scattering distributes particles over 
several canonical orbitals located around the Fermi level. 
As compared to the implementation of the EDF based on 
a Slater determinant, this might lead to the promotion 
of particles from low/high angular-momentum states to 
high/low angular momentum orbitals [t^ . Depending on 
the situation, this will favor or inhibit the formation of 
halos. As opposed to the anti-halo effect discussed above, 
the way this process impacts halos depends on the sys- 
tem and on the particular distribution of orbitals around 
the Fermi energy at the drip-line. 



III. EXISTING INVESTIGATIONS AND 
ANALYSIS METHODS 



C. Role of pairing correlations 

Theoretical investigations of nuclei far from stability, 
either within non-relativistic [63 . Issl . [66| or relativis- 
tic [13, m, [69|] EDF frameworks, have pointed out the 
importance of pairing correlations. This makes the im- 
plementation of the SR-EDF method in terms of a quasi- 
particle vacuum more successful in practice than the one 
based on a reference Slater determinant. 

The explicit treatment of pairing correlations through 
dependencies of the nuclear EDF on the anomalous den- 
sity changes qualitatively the density profile in loosely 
bound systems. By studying the asymptotic form of 
the quasiparticle wave- functions solution of Eq. ([5]), 
it is easy to show that the decay constant kq at 
play is now kq — \J ~2meQlW- , where \tQ\= Eq — \ and 
Eq = mini/[£',y] is the lowest quasiparticle energy solu- 
tion of Eq. ([5|). Considering the most extreme case 
of a canonical state lying at the Fermi level at the 
drip-line (eg ~ A « 0), one sees that |eo| ~ i^o ~ > 0. 
Therefore, everything else being equal, paired densities 
decrease faster than unpaired ones at long distances. 
Because the decay constant does not go to zero as 
eo « A ~ 0, the second moment of the density cannot 



Halo properties of medium-mass drip-line nuclei have 
been studied for various isotopic chains u sing relativistic 
or non-relativistic EDF methods [Ei, III, [la III, [t^, S 
[t^, [13, HH . Owing to the discussion provided above, the 
evolution of the r.m.s radii along isotopic chains is often 
used to characterize halos in a qualitative manner. One 
needs however more quantitative characterizations of the 
halo itself. For example, the concept of giant halo was 
recently introduced on the basis of summing up the occu- 
pations of low-lying orbitals with large r.m.s. radii [tsI ]. 
Such halo structures, supposedly composed of six to eight 
neutrons, have been characterized through relativistic 
and non-relativistic methods [tI, [t^, [t^ Ho, [81| , mainly 
for Zr and Ca isotopes, and were related to the presence 
of £ — 1 states close to the Fermi level at the drip-line. 
Finding giant halos in medium-mass nuclei is intuitively 
surprising. Indeed, spatially decorrelated neutrons seem 
less likely to appear as the mass of the system increases 
and their behavior tends to become more collective. We 
will come back to this point. 

The present section is devoted to discussing observ- 
ables and analysis tools that are usually used to identify 
and quantify halo signatures in nuclear systems. The 
purpose is to introduce generic features which turn out 
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to be useful later on and, above all to demonstrate the 
limitations of existing analysis tools. 

Chromium and tin isotopic chains are chosen as test- 
ing cases throughout this work. Calculations are per- 
formed using the non-relativistic HFB spherical code 
HFBRAD [13 • In HFBRAD, the space is discretized within 
a sphere using vanishing boundary conditions for the 
wave functions (Dirichlet conditions). Convergence of the 
calculations as a function of numerical parameters has 
been checked for all results presented here. The Skyrme 
SLy4 functional [H, HI] is employed in the particle- 
hole channel. The particle-particle effective vertex is a 
density-dependent delta interaction (DDDI) correspond- 
ing to a "mixed- type" pairing. Its density-dependent 
form factor is a compromise between a pairing which is 
constant over the nucleus volume ("volume- type"), and 
one which is peaked at the nucleus surface ("surface- 
type") Hi, Hi, 113, m. Hi]. To avoid the ultra-violet 
divergence associated with the local nature of the pair- 
ing functional, a phenomenological regularization scheme 
corresponding to a smooth cutoff at 60 MeV in the single- 
particle equivalent spectrum is used [65j . Such a pairing 
scheme is referred to as REG-M. 

The HFB problem is solved self-consistently. Thus, 
the shape of the central potential cannot be manually 
adjusted to reduce the binding energy of weakly-bound 
orbitals and halo candidates can only be identified a pos- 
teriori. 



A. First characterizations 

1. Chromium isotopes 

Among all medium-mass nuclei predicted to be 
spherical [o^, [Mli Chromium isotopes {Z = 24) located 
at the neutron drip-line are good halo candidates. In 
Fig. [H neutron canonical energies e" in the vicinity 
of the positive energy threshold are plotted along the 
Chromimn chain, ^''Cr being the predicted drip-line 
nucleus. The presence of low-lying 3si/2 and 2^5/2 
orbitals at the drip-line provides ideal conditions for the 
formation of halo structures. 



As discussed in Sec. Ill Bl the abnormal extension of 
the one-body neutron density is usually characterized 
through the evolution of the neutron r.m.s. radius as 
one approaches the drip-line, as presented in Fig. [31 A 
significant kink in the neutron r.m.s. radius is seen at the 

= 50 shell closure. Such a kink is usually interpreted 
as a signature of the emergence of a neutron halo [H, [t^ . 
However, this could equally be due to a simple shell effect. 
Indeed, as the A^ = 50 gap is crossed, the two-neutron 
separation energy S2n drops, as seen in Fig. [H As a re- 
sult, the decay constant kq of the one-body density is 
largely reduced. However, a genuine halo phenomenon 
relates more specifically to the presence of nucleons which 
are spatially decorrelated from a core. Even though the 




Figure 1: (Color Online) Neutron canonical energies e" along 
the Cr isotopic chain, obtained through spherical HFB cal- 
culations with the {SLy4-|-REG-M} functional. Conventions 
used in all the figures of the article are given in Fig. [2l 
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Figure 2: (Color Online) Conventions used in all the figures 
for the labeling of individual states and of the chemical po- 
tential. 



case of drip-line Cr isotopes seems favorable, as the S2n 
drops to almost zero at A^ = 50, the occurrence of a halo 
cannot be thoroughly addressed by only looking at the 
evolution of the neutron r.m.s. radius. 



2. Tin isotopes 

Sn isotopes {Z — 50) are considered as a milestone 
for EDF methods and are rather easy to produce in ra- 
dioactive beam experiments because of their magic pro- 
ton number. In particular, the fact that it is a long 
isotopic chain is convenient for systematic studies. At 
the neutron drip-line, which corresponds to ^^^Sn for the 
{SLy4-|-REG-M} parameter set, the least-bound orbitals 
are mostly odd-parity states. Among them, 8^3/2 and 
3pi/2 states might contribute significantly to the forma- 
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Figure 3: (Color Online) Same as Fig. [1] for proton (■), 
neutron (♦) and charge (A) r.m.s. radii. Experimental values 
for charge r.m.s. radii are indicated when available (t), along 
with experimental error bars [o^ . 
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Figure 4; (Color Online) Same as Fig.[T]for two-neutron sep- 
aration energies S2n (♦). Experimental values are indicated 
when available (^ when both masses are known, T when 
at least one comes from mass extrapolation), along with ex- 
perimental error bars. 



tion of a halo (Fig. [5]). 

However, whereas those I — \ states are relatively well 
bound, the least bound orbital is the l«i3/2 = 6) in- 
truder state which is strongly affected by the confining 
centrifugal barrier. Nevertheless, the neutron r.m.s. ra- 
dius (Fig.[6|) exhibits a weak kink at the iV = 82 shell clo- 
sure, which has been interpreted as a halo signature [5^. 

As pointed out previously, an analysis based only on 
r.m.s. radii is incomplete and can even be misleading. 
Indeed, although the shell effect at the = 82 magic 
number generates a sudden decrease of the S2n, the lat- 
ter does not drop to zero, as seen in Fig. [71 A direct 
connection between the kink of the r.m.s. radius and the 
formation of a neutron halo is thus dubious. This point 
will be further discussed below. 
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Figure 5: (Color Online) Same as Fig. [T]for Sn isotopes. 
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Figure 6: (Color Online) Same as Fig. [3]for Sn isotopes. 



In any case, the analysis based on neutron r.m.s. radii 
is useful but insufficient to characterize halo in a man- 
ner that allows the extraction of information useful to 
nuclear structure and theoretical models. The charac- 
terization of halos through the definition of the neutron 
matter thickness and the one-neutron region thickness 
is possible but remains arbitrary and correlated to 
a one-neutron halo hypothesis. Another way is to ex- 
tract so-called "halo factors" from the individual spec- 
trum through antiproton annihilation probing the nu- 
clear density extension JMJ. However, such tools do 
not allow the extraction of quantitative properties, such 
as the actual number of nucleons participating in the 
halo. They also define the halo as the region where the 
neutron density dominates the proton one, which is an 
admixture of the neutron skin and the (potential) halo. 
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Figure 7: (Color Online) Same as Fig. |4] for Sn isotopes. 



B. The Helm model 

1. Introduction 

The Helm model has recently been exploited to rem- 
edy to the lack of quantitative measure of halo existence 
and properties 59). Originally, the purpose of the Helm 
model Refs. (95j, i96|, i97j] was to fit experimental charge 
densities, using a few-parameter anzatz, in view of ana- 
lyzing electron scattering data. The normalized nuclear 
charge density is approximated by the convolution of a 
sharp-sphere density of radius Rq defining the nuclear 
extension and of a gaussian of width cr describing the 
surface thickness. The r.m.s. radius of the Helm density 
solely depends on Rq and cr and reads as 

This model has been used to study neutron skins and 
halos in medium-mass nuclei close to the neutron drip- 
line (soj . Proton and neutron densities were defined as 
a superposition of a core density pl^^^ plus a tail density 
/^taii describing, when appropriate, the halo. The idea 
was to reproduce the core part p^^^^ using the Helm an- 
zatz normalized to the nucleon number N'^ (A^^ — N 
or Z). Thus, the two free parameters {Rq,(j'^) were ad- 
justed on the high momentum part of the realistic form 
factor 



= 47r y p'?(r) jn{k r) dr , (14) 

where p'^(r) is the density coming out of the many-body 
calculations. It was suggested in Ref. [5§| to evaluate 
(i) Rq through the first zero kf of the realistic form fac- 
tor: Rq — zl/kf, where zl is the first zero of the Bessel 
function ji (zl « 4.49341), and (ii) cr* by comparing the 



model and realistic form factors at their first extremum 
kfj (a minimum in the present case). Then, the fol- 
lowing radii are defined (i) i?geom('?) ~ ^5/3 i?r.m. 3.(9) 
(geometric radius) for realistic densities, and (ii) 

i?Hcim(g) = ^/5/3R^„,,,iq) = ^jRf + 5ai^ (Helm ra- 
dius) for model densities. 

Adjusting the Helm parameters to the high momentum 
part of the realistic form factor was meant to make the 
fitting procedure as independent of the asymptotic tail 
of p'^{r) as possible. Constructed in this way, i?Hoim(«) 
should not incorporate the growth of i?geom("-) when the 
neutron separation energy drops to zero and the spa- 
tial extension of weakly-bound neutrons increases dra- 
matically. In addition, it was checked that the difference 
between i?Hcim(p) and -Rgoom(p) was negligible near the 
neutron drip-line. From these observations, the neutron 
skin and neutron halo contributions to the geometric ra- 
dius were defined as^^^ 



Ai?skin — R}ichn{n) — i?Hclm (p) , 

(15) 

^^halo — -^gcom 



2. Limitations of the Helm model 

Proton and neutron Helm radii are compared to ge- 
ometric ones on Fig. [5] for chromium and tin isotopes. 
The behavior of i?geom(9) and i?Hcim(9) for Sn isotopes 
is the same as in Ref. [H^*^''-'. For both isotopic chains, 
the sudden increase of the neutron geometric radius be- 
yond the last neutron shell closure might be interpreted 
as a signature of a halo formation. However, Ai?haio is 
non-zero along the entire Cr isotopic chain, even on the 
proton-rich side. The latter result is problematic as neu- 
tron halos can only be expected to exist at the neutron 
drip-line. 

Such non-zero values for A_Rhaio can be understood 
as a direct consequence of the gaussian folding in the 
definition of the Helm density. The asymptotic decay 
of the Helm density is roughly quadratic in logarithmic 
scale, instead of being linear [6J, [H, . To illustrate 
this point, Fig.[5]displays the realistic and Helm densities 
of ^"'Cr (in the valley of stability) and ^''Cr (drip- line 
nucleus). The difference in the asymptotic behaviors is 
obvious. In particular, the Helm densities are unable 
to reproduce the correct long-range part of the non-halo 



Similar definitions could be applied to nuclei close to the proton 
drip-line, where a proton halo is expected instead of a neutron 
one. 

* Results diff'er slightly from Ref. [59j because of the different pair- 
ing functional and regularization scheme used, as well as the 
larger number of j-shells taken into account in the present cal- 
culations. 
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Figure 8: (Color Online) Geometric and Helm radii for Cr 
and Sn isotopes calculated in the spherical HFB approach 
with the {SLy4+REG-M} functional. 



proton density, or the neutron density of nuclei in the 
valley of stability. 

Such features lead to unsafe predictions for the halo 
parameter Ai?haio as the neutron skin and the potential 
halo are not properly separated. Such problems, as well 
as a lack of flexibility to account for finer details of the 
nuclear density had already been pointed out [99] . 

One might thus question the fitting procedure intro- 
duced in Ref. [5§|. The method naturally requires Rq 
and cr' to be adjusted on the form factor at sufficiently 
large k so that the Helm density relates to the core part 
of the density only. Of course, some flexibility remains, 
e.g. one could use the second zero fc| of F'(fc) to adjust 
Rq. Following such arguments, four slightly different ad- 
justment procedures Ai, i = 1,4, all consistent with the 
general requirement exposed above, have been tested to 
check the stability of the Helm model 



Ai : (i) 
A2 : (i) 

A, : (i) F^(fc|) 



F{kf) {ii) Ff,{kl,)^F^{k[ 



M) ' 




4.0 8.0 

r [fm] 



Figure 9: (Color Online) Realistic and Helm densities of ^*Cr 
and '^''Cr. 



A, : (i) F^^ikl) ^ F'^ikl) (ii) Ffj'iOAkD^Fi'iOAkl). 

Fig. [TO] shows the halo parameter Ai?haio obtained for 
Cr isotopes using protocols Ai to A4. Note that proto- 
col Ai is the one proposed in Ref. [1^ and used earlier 
whereas the weight of the long distance part of the real- 
istic density is more important in protocol A^. Although 
the general pattern remains unchanged, the halo param- 
eter significantly depends on the fitting procedure used 
to determine (i?Q,(7'^). Because of the wrong asymptotic 
behavior of the Helm density discussed above, one can- 
not make Ai?haio to be zero for magic and proton-rich 
nuclei (see protocol ^4), keeping unchanged its values 
for halo candidates at the neutron drip-line*^^^ Such a 
fine tuning of the fitting procedure that would make use 
of an a priori knowledge of non-halo nuclei is impractical 
and unsatisfactory. 

As a next step, we tried to use other trial densities to 
improve on the standard Helm model. A key feature is to 



FHkf) (iiX(fc?)=F'^'(fc?), 
FHkl) {:n)Ff,'{kl)=F<i'{kl). 



Helm densities obtained with the Ai protocol still do not match 
the realistic ones, even for protons. 
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Figure 10: (Color Online) Halo parameter Ai?haio for 
chromium isotopes using different fitting procedures for the 
Helm parameters {Rq,(j'') (see text). 



obtain an analytical expression of the associated model 
form factor in order to adjust easily its free parameters. 
We could not find any form leading to both an analytical 
expression of F'^ and good asymptotic, with only two 
free parameters'^''. 

Although the Helm model looked promising at first, 
we have shown the versatility of its predictions. The 
inability of the model to describe the correct asymptotic 
of the nuclear density in the valley of stability, as well 
as the too large freedom in the fitting procedure, limit 
very much its predictive power. Therefore a more robust 
analysis method is needed to characterize medium-mass 
halo nuclei. 



IV. NEW CRITERION FOR A QUANTITATIVE 
ANALYSIS OF HALO SYSTEMS 

Although deceiving, the previous attempts have under- 
lined the following point: a useful method to study halos 
must be able to characterize a spatially decorrelated com- 
ponent in the nucleon density in a model-independent 
fashion. We propose in the following a method which al- 
lows the identification of such a contribution to the inter- 
nal one-body density. Our starting point is a thorough 
analysis of medium-range and large-distance properties 
of the one-body internal density in Sec. IIV Al Based on 
such an analysis, new quantitative criteria to identify and 
characterize halos are defined in Sec. IIVBI We already 
outline at this point that the analysis and the associated 
criteria are applicable to any finite many-fermion system. 



Using model densities depending on three parameters would 
make the Helm model even more dependent on the fitting pro- 
cedure. 



as long as the inter-fermion interaction is negligible be- 
yond a certain relative distance. As done throughout the 
article, atomic nuclei are used as typical examples in the 
present section. 

A. Properties of the one-body density 

1. Definitions and notations 

Let US start from the non-relativistic iV-body 
Hamiltonian'^) 

N 2 ^ 

^"-El^+E^-^.)' (16) 

i—l — 1 

i<j 

where pi is the single-particle momentum, = l^'i — | 
and V denotes the vacuum nucleon-nucleon interaction. 
The nuclear Hamiltonian H'^ is invariant under trans- 
lation and can be written as a sum of a center-of-mass 
part iJ^jjj and an internal part H^^. Thus, eigenstates 
of , denoted by '^fgi'^'i ■ ■ - ^n), can be factorized into 
the center-of-mass part (plane wave) times the internal 
wave function 

^f^^in ■ ■ ■ tn) = e'^-^" I'f (a • • • 6v-i) , (17) 

where K is the total momentum and the center-of- 
mass position 



1=1 

The word internal relates to the fact that the wave func- 
tion can be expressed in terms of relative coordinates 
only, such as the {N — 1) independent Jacobi variables 



ii = n+i--E^j' (1^) 

and is associated with the internal energy . A con- 
sequence is that is invariant under translation of the 
system in the laboratory frame. 

The ground-state internal wave function $^ can be 
expanded in terms of the complete orthonormal set of 



The Coulomb interaction is omitted here, as the focus is on neu- 
tron halos. The spin degrees of freedom are also not explicitly 
included as their introduction would not change the final results. 
Finally, the Hamiltonian is restricted to a two-body interaction. 
The conclusions would not change either with the introduction 
of the three-body force. 
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internal (TV — l)-body wave functions which 
are e igen s tate s of the (N — l)-body internal Hamilto- 
nian [lo3, [IqI EM, flO^ 



with = 2mi.cdeiy/^^, and e^, = (i?^ — E^~^^ is 
minus the one-nucleon separation energy to reach 
Solutions of the free Shrodinger equation take the form 



(20) 



such that 



<(rl . . . r-^) = ^ <-i(fi . . . rr,-i) 

X ip^{fN - Rn^i) ■ (21) 



The states <&^~^ are ordered by increasing ener- 
gies, 1/ — corresponding to the ground state of the 
{N — l)-body system. The norm of the overlap func- 
tions (£j/(r^)provides the so-called spectroscopic fac- 
tors [lM . fl05| 



(22) 



Finally, the relevant object to be defined for self-bound 
syst ems is the internal one-body density matrix [102i . il06l 



Pli]ir,r') = ^ (/3*(r ')¥>,. (r ) , 



(23) 



which is c ompletely determined by the overlap func- 
tions [l06l |. The actual internal one-body density 
P[i] (^) = P[i] ('') is extracted as the local part of the in- 
ternal density matrix 



P[i](r)^5^|^.(f)p=5]?^|^.(r)p 



An 



(24) 



where the energy degeneracy associated with the orbital 
momentum has been resolved through the summation 
over the spherical harmonics. 



2. Long-distance behavior and ordering of the ipv{f) 

For large distance, i.e. r > R, the nuclear interaction 
vanishes and the asymptotic radial part (fu^^^ of the over- 
lap function is solution of the free Schrodinger equation 
with a reduced mass mrcd = m{N — 1)/N 



(P ^2 d C(C + 1) 
dr^ r dr 



^^M-O, (25) 



* In the following, the radial part of a wave function /(r) is noted 



(26) 



As a result, the internal one-body density behaves at long 
distances as^^-' 



(27) 



For very large arguments, the squared modu lus of a 



Hankel function behaves as e 



-2Kir 



Thus 



the = component dominates and provides the usual 
asymptotic behavior [H [H, [oll^i") 



(Kor)2 



(28) 



The asymptotic form of the Hankel function is inde- 
pendent of the angular momentum which explains why 
high-order moments (r") of the density diverge when 
high-^ states are loosely bound, as discussed in Sec. IIIBI 
Thus, the contributions of the overlap functions to pj^j at 
very large distances are ordered according to their asso- 
ciated separation energies |e^|, independently of ly. Cor- 
rections to this ordering at smaller distances come from 
(i) the ^-dependence of the Hankel functions due to the 
centrifugal barrier, which favors low angular momentum 
states, and (ii) the {21 -I- 1) degeneracy factor which fa- 
vors high angular momentum states. In any case, for ex- 
tremely large distances the least bound component will 
always prevail, although this may happen beyond simu- 
lation reach. 

To characterize the net effect of corrections (i) and (ii) 
on the relative positioning of overlap functions at long 
distances, the contributions {2iy -f l)|(^j,(r)p, for a fixed 
energy but different angular momenta, are compared in 
Fig. [11] for the solutions of a simple finite spherical well. 
Outside the well, Hankel functions are exact solutions of 
the problem. The potential depth is adjusted to obtain 
identical eigenenergies for all 1^. Although the (2^^ -I- 1) 
factor reduces the gap between s and p components, the 
effect of the centrifugal barrier is always the strongest 



Rigorously, this is true only if the convergence of the overlap 
functions to their asymptotic regime is uniform in the mathemat- 
ical sense, i.e. if they reach the asymptotic regime at a common 
distance R [lQ6]. This is not actually proven in nuclear physics, 
but it has been sh own to be true for the electron charge density 
in atomic physics [lOSL Il09j . 

Note that the asymptotic of pP and p" are different because 
of the charge factor (Hankel functions for neutrons, Whittaker 
functions for protons). 
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at large r, where states are clearly ordered according to 
favoring low angular momenta. In any case, the sep- 
aration energy remains the leading factor as far as the 
ordering of overlap functions at long distances is con- 
cerned. 




r [fm] 

Figure 11: (Color Online) Squared components of the solu- 
tions of a finite spherical well of fixed radius a = 4 fm, multi- 
plied by the degeneracy factor (2^^ -I- 1), for various angular 
momenta and fixed energy — —100 keV. The first state 
for each £^ (nodeless component corresponding to a primary 
quantum number equal to zero) is represented. 



more bound. Of course, the centrifugal barrier influences 
the position of such crossings but not their occurrence be- 
cause of the robustness of the (very) asymptotic ordering 
pattern discussed in the previous section. 

Let us now incorporate the role of spectroscopic fac- 
tors. In practice, Si, is known to increase with the ex- 
citation energy of the corresponding eigenstate of the 
(TV — l)-body system. Thus, the norm of ipo is smaller 
than those of the excited components ip^, which mechan- 
ically ensures the existence of the crossings discussed pre- 
viously. A similar reasoning holds when going from ipo 
to (fi etc. 

One should finally pay attention to the number of 
nodes of the overlap function (p,^. This feature actu- 
ally favors low angular momentum states as far as the 
asymptotic positioning is concerned. If two components 
have the same energy but different angular momenta, 
the one with the lowest £ will have a greater number of 
nodes. This will reduce the amplitude of the wave func- 
tion in the nuclear interior. That is, the weight of the 
asymptotic tail is increased, which favors its dominance 
at long distance. However, this effect is expected to have 
a small impact in comparison with the other corrections 
discussed above. As a result, the crossing pattern be- 
tween the components of the density is not jeopardized 
by the existence of nodes in the overlap functions. 



B. Halo characterization 



3. Crossing pattern in (r) 



1. Defimtion 



The (model-independent) ordering at long distances of 
individual components entering pj^j has interesting con- 
sequences on the properties of the density as a whole. As 
discussed below, this ordering induces a typical crossing 
pattern between the individual components which will 
eventually be used to characterize halo nuclei. 

Introducing normalized overlap functions ip^^r), 
Eq. pi)) becomes 



(29) 



Let us take all spectroscopic factors equal to one for 
now. The v = ^ component, corresponding to the small- 
est separation energy, dominates at large distances. Be- 
cause of continuity and normalization conditions, this im- 
plies that ?/'o(?') has to cross all the other overlap func- 
tions as r goes inward from +cx3 to zero. The position 
at which -00 crosses each t/j^ depends on the difference of 
their separation energies and on their angular momenta. 
In particular, there will exist a crossing between |'!/jo(T')p 
and the remaining density {p^^^if) - Co(r)] . The same is 
true about |i/;i(r)p: it must cross the remaining density 
[p[^l(r) — Co^r) — Ci(r)] ... As a result, any given indi- 
vidual component must cross the sum of those that are 



The discussion of Sec. IIV A 31 demonstrates how indi- 
vidual contributions to the one-body density (i) are po- 
sitioned with respect to each other (ii) display a typical 
crossing pattern. Such features are now used to charac- 
terize halo systems. 

As pointed out earlier, one general and model- 
independent definition of a halo relates to the existence 
of nucleons which are spatially decorrelated from oth- 
ers, constituting the core. This can only be achieved if 
some contributions to the internal density exhibit very 
long tails. Most importantly, the delocalization from the 
core requires the latter to exist and to remain well lo- 
calized. To achieve such a spatial decorrelation between 
a core and a tail part, it is necessary to have a crossing 
between two well-identified groups of orbitals with signif- 
icantly different asymptotic slopes. This translates into a 
sharp crossing between those two groups of orbitals and 
thus to a pronounced curvature in the density. Note that 
this explains the empirical observation that the first log- 
arithmic derivative of th e de nsity invariably displays a 
minimum at some radius llll| . How much this feature is 



pronounced or not is key and will be used in the follow- 
ing to design model-independent criteria to characterize 
halo systems. 

A pronounced crossing is illustrated in Fig. [12] for a 
simple model where the halo is due to a single orbital. 
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Of course, more complex situations have to be considered 
where multiple states contribute to the core and the halo. 
Indeed, the presence of collective motions in medium- 
mass systems implies that one hardly expects a single 
orbital to be well separated from the others. 




10 15 

r [fm] 



Figure 12: (Color Online) " Core+tail" simplified model. The 
total density is the superposition of a well-bound component 
and a loosely-bound one. A semi-phenomenological density 
(see Appendix[Bj is used for the core density, whereas the halo 
part is the realistic 3i/2 state of *'^Cr obtained from spherical 
HFB calculations with the {SLy4-|-REG-M} functional. 



2. Relevant energy scales 

The need for the existence of two groups of orbitals 
characterized by significantly different asymptotic slopes 
provides critical conditions for the appearance of a halo: 
(i) the least bound component ipo must have a very small 
separation energy to extend far out, (ii) several compo- 
nents ipi^ip2 ■ ■ ■ ifm may contribute significantly to the 
density tail if, and only if, they all have separation ener- 
gies of the same order as that of (fa, (iii) for this tail 
to be spatially decorrelated from the rest of the den- 
sity (the "core"), the components with v > Vm have to 
be much more localized than those with u < Vm- This 
third condition is fulfilled when the crossing between the 
m}^ and (m -I- 1)*'' components in the density is sharp, 
which corresponds to significantly different decay con- 
stants Km ^ tm+1 at the crossing point. 

The later situation translates eventually into spe- 
cific patterns in the excitation energy spectrum of the 
(TV — l)-body system. It suggests that a halo appears 
when (i) the one- neutron separation energy S'„ = |eo| is 
close to zero, (ii) a bunch of low energy states in the 
(iV — l)-body system have separation energies close 
to zero, and (iii) a significant gap in the spectrum of the 
[N — l)-body system exists, which separates the latter 
bunch of states from higher excitations. 

A similar discussion was given in the context of de- 
signing an effective field theory (EFT) for weakly-bound 



nuclei [112l |. where two energy scales {E,E') were found 
to be relevant: (i) the nucleon separation energy E — Sn 
which drives the asymptotic behavior of the one-body 
density, and (ii) the core excitation energy E' — |em+i| 
which needs to be such as E' ^ E, in order for the tail 
orbitals to be well decorrelated from the remaining core. 
The additional energy scale that we presently identify 
is the energy spread AE' of the low-lying states in the 
{N ~ l)-body system, which becomes relevant when more 
than one component is involved in the halo. The cor- 
responding picture is displayed in the bottom panel of 
Fig. [T2] and is also translated in terms of canonical ener- 
gies Ci in the upper panel of the same figure. 
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(a) Canonical neutron energy spectrum 




(b) Separation energy spectrum |e^| for the (A'' — l)-body system. 



Figure 13: (Color Online) Schematic display of the energy 
scales relevant for the appearance of halos (right-hand side). 
The realistic spectra obtained through HFB calculations of 
the four last bound chromium isotopes are shown on the left- 
hand sides. 

More quantitatively, the ideal situation for the forma- 
tion of a halo is obtained for (i) a very small separation 
energy, in orders of a few hundred keVs, the empirical 
value of 2 MeV/A^/^ from Refs. [H, [13 giving a good 
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approximation of expected values, (ii) a narrow bunch 
of low- lying states, whose spread AE should not exceed 
about one MeV, and (iii) a large gap E' with the remain- 
ing states, at least four or five times the separation energy 
E. Those are only indicative values, knowing that there 
is no sharp limit between halo and non-halo domains. 



3. Halo region 

As discussed in the previous section, a halo can be 
identified through a pronounced ankle in the density, due 
to the sharp crossing between the aggregated low-lying 
components and the upper- lying ones. Such a large cur- 
vature translates into a peak in the second derivative of 
the (base-10) logarithmic profile (logj^g) of the one-body 
density, as seen in Fig. [14] for a schematic calculation. 
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the halo, if it exists, dominates. Therefore, we define the 
spatially decorrelated region as the region beyond the ra- 
dius ro where the core density is one order of magnitude 
smaller than the halo one. In practice, the previous defi- 
nition poses two problems. First, in realistic calculations, 
one only accesses the total density. Second, the choice of 
one order of magnitude is somewhat arbitrary. 

Extensive simulations have been performed to charac- 
terize ro unambiguously, using either one or several con- 
tributions to the halo density, and covering large energy 
ranges for E, E' and AE. More details on the method 
used to find the best approximation to rg, as well as the 
corresponding theoretical uncertainty, are given in Ap- 
pendix |B] Given rmax , which can be extracted from the 
total density, it has been found that tq can be reliably 
defined through 



' max I 

9^ logio P{r) 



Qj.2 



r = ro 



^2dHog,,p{r) 
5 



(30) 

as exemplified in Fig. [151 Also, theoretical uncertainties 
on the determination of ro are introduced, such that 



0.35 <^^|^^#^< 0.50, 
where ' denotes a compact notation for d/dr. 



(31) 



Once validated by simulations, the method to isolate 
the halo region only relies on the density as an input, and 
does not require an a priori separation of the one-body 
density into core and halo parts. Finally, one may note 
that our definition of the halo region does not a priori 
exclude contributions from individual components with 
angular momenta greater than one. 



4. Halo criteria 

We now introduce several criteria to characterize the 
halo in a quantitative way, by applying the previous anal- 
ysis to the neutron one-body density'^^^^ First, the av- 
erage number of nucleons in the halo region can be ex- 
tracted through 



Figure 14: (Color Online) Ankle in the (base-10) log-density 
due to the presence of a low-lying state well separated from 
the remaining ones: log-density (bottom panel), first (middle 
panel) and second (top panel) log-derivatives. The conven- 
tions are the same as in Fig. 1121 



halo 



r+00 

in / p" (r) r^ dr 



(32) 



At the radius r = r^nax corresponding to the maximum 
of that peak, core and tail contributions cross; i.e. they 
contribute equally to the total density. At larger radii. 



For neutron-rich medium-mass nuclei, protons are well confined 
in the nuclear interior, thus do not participate in the long-range 
part of the total density p. The two densities p and p" can be 
used regardless to evaluate A'^ij^io E'nd i5_Rijaio- 
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Figure 15: (Color Online) Definition of ro through the second 
derivative of the log-density, using the same model density as 
in Fig. [ill '"0 is represented by the central vertical line. The 
shaded area corresponds to the tolerance margin on ro (see 
text). 



The quantity (5i?haio is similar to Ai?haio defined within 
the Helm model (Eq. p3)) ). However, the former does 
not rely on any a priori decomposition of the density into 
core and halo components. That is of critical importance. 
Extensions to all radial moments of the density can be 
envisioned*^^^^. 

The quantities A'haio and SRhsdo are of course corre- 
lated, but they do not carry exactly the same informa- 
tion. Note that tolerance margins on tq from Eq. (j3ip 
propagate into theoretical uncertainties on A'^iiaio and 

In the case of stable/non-halo nuclei, both quantities 
will be extremely small. There is still a slight curvature 
in the density profile that provides a radius rg but the 
computed criteria will be consistent with zero. In the 
particular case of magic neutron number, the curvature 
becomes particularly weak and translates into a broad 
peak in the second log-derivative. As a result, the radius 
To value is large and defines a region where the density is 
particularly low. This is illustrated by Fig. 1171 where rg 
is plotted for chromium isotopes as a function of A. The 
maximum of ro is attained for the magic shell N = 50. 
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Figure 16: (Color Online) Consequences of the definition of 
ro (vertical lines for the values of ro and the tolerance margin 
- see text) in the same model as in Fig. 1141 The halo density 
dominates the core part by around one order of magnitude. 



An important information is the effect of the halo re- 
gion on the radial moments of the density. By definition, 
the contribution of the core to any moment (r") is neg- 
ligible for r > ro. It has been checked in the case of the 
r.m.s. radius, and is all the more true as n increases. 
Thus, one can evaluate the effect of the decorrelated re- 
gion on the nuclear extension through 



Figure 17: (Color Online) Evolution of ro along the Cr iso- 
topic chain, obtained through spherical HFB calculations with 
the {SLy4-|-REG-M} functional. 

Finally, further characterization of the halo can be 
achieved by looking at the individual contributions of 
each overlap function 
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(33) 



Numerical issues appear when going to high-order moments. In- 
deed, (r") is more and more sensitive to the upper limit of inte- 
gration as n increases. Thus, the result may significantly depend 
on the box size used to discretize the continuum or on the size 
of the basis used to expand quasiparticle wave-functions in HFB 
calculations. 
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Miaio,i/ provides a decomposition of the halo in terms of 
single-particle-hke states. Note that the inner part of an 
overlap function, i.e. for r < tq, does not contribute to 
halo observables. 

By analogy with the criterion used for light halo sys- 
tems, the probability of each individual overlap function 
to be in the r > tq region can be defined through 

/ + 00 
|(^.(r)pr2dr 

t^^^-^ ■ (35) 

/ |(^,(r)|2r2dr 

V. APPLICATION TO EDF CALCULATIONS 

We apply the analysis method introduced in Sec. IIVI 
to results obtained from self-consistent HFB calculations 
of chromium and tin isotopes. In Sec. IIVI the energies el 
that characterize internal overlap functions denote exact 
nucleon separation energies. No approximation to the 
nuclear many-body problem was involved in the analysis 
conducted in Sec. IIVI The patterns of the internal one- 
body density thus extracted are fully general and model- 
independent. 

In practice of course, one uses an approximate treat- 
ment of the quantum many-body problem. This raises 
critical questions in the case of EDF calculations as dis- 
cussed in Appendix IA2I Indeed, the one-body density 
at play in single-reference EDF calculations is an intrin- 
sic density rather than the internal density, i.e. it is the 
laboratory density computed from a symmetry breaking 
state. As is customary in EDF methods though, one uses 
such an intrinsic density to approximate the internal den- 
sity; e.g. when analyzing electron scattering data. Of 
course, such an identification is not rigourously justified 
and formulations of EDF methods directly in ter ms o f 
the internal density are currently being considered |ll3l | . 
Still, the asymptotic part of the lower component VJ{r) 
of the HFB quasiparticle wave-function satisfies the free 
Schrodinger equation [sll (Eq.[25]), just as the true inter- 
nal overlap function ^Pi,{r) does. Considering in addition 
that the intrinsic HFB one-body density reads as 

^.(^)^^M±l|v;.(,)|2^ (36) 

one realizes that the analysis performed in Sec. IIV| in- 
cluding the existence of the crossing pattern, applies di- 
rectly to it("). 



A. Implementation of the criteria 

In the code HFBRAD, the HFB problem is solved in a 
spherical box up to a distance Rhox from the center of 
the nucleus on a radial mesh of step size Ar = 0.25 fm. 
For Rbox = 40 fm, the mesh has 160 points in the radial 
direction, for both the individual wave-functions and the 
densities. To obtain a satisfactory precision, the second 
order log-derivat ive i s computed using a five-points dif- 
ference formula [ll0| |. The precision of the formula is 
the same as the intrinsic precision of the Numerov algo- 
rithm used for the integration of s econ d-order differential 
equations (which is O (Ar^)) [83 . lll4l |. Approximate po- 
sitions of the maximum of the second order log-derivative 
of p" (r) and of tq are first determined with a simple com- 
parison algorithm. To increase the precision, a 11-points 
polynomial spline approximation to the density and its 
second log-derivative around the two points of interest 
is performed. Because the functions involved are regular 
enough, a spline approximation provides the radii rmax 
and To with a good precision, as they are obtained us- 
ing a dichotomy procedure up to a (arbitrary) precision 
of 10^^. Finally, the integrations necessary to compute 
A^haio and 6Rhaio are performed with six-points Gaussian 
integration. 

In the definition of (5i?haio, the core contribution to the 
total r.m.s. radius is approximated as the root-mean- 
square radius of the density distribution truncated to its 
r < tq component. To check the influence of this cut, the 
core density was extrapolated beyond the point where the 
second order log-derivative crosses zero*-^''-' using Eq. ((28)) 
and enforcing continuity of p" and p"'. No difference was 
seen for 5i?haio- 

The individual contributions iVhaio,i, as well as the in- 
dividual probabilities Pi, are evaluated in the canonical 
basis. Equivalently, iVhaio.i/ and can be calculated in 
the quasiparticle basis. Quasiparticle states are the best 
approximation to the overlap functions, but canonical 
and quasiparticle basis really constitute two equivalent 
pictures. Indeed, each canonical state is, roughly speak- 
ing, split into quasiparticle solutions of similar energies. 
A summation over quasiparticles having the same quan- 
tum numbers in an appropriate energy window would 
recover the single-particle canonical approximation. The 
latter is preferred here, as it is more intuitive to work in 
the natural basis. 



B. Cr isotopes 

According to the analysis of Sec. IIVB 21 drip-line 
chromium isotopes appear to be ideal halo candidates. 



The method was developed in Sec. llVl for the exact internal den- 
sity in order to demonstrate its generality and to eventually ap- 
ply it to the results of other many-body methods dealing with a This is the point where the halo contribution effect becomes sig- 
variety of finite many-fermion systems |58ll . nificant. 
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The separation energy spectrum |e"| = — A" to the 
states in the {N — l)-body system is shown in Fig. [THl 
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Figure 18: (Color Online) Neutron separation energies 
|e"| = E" — A" along the Cr isotopic chain, obtained through 
spherical HFB calculations with the {SLy4+REG-M} fiinc- 
tional. Only relevant quasiparticle energies (A^" > 0.01) are 
displayed. Conventions for labeling individual states are 
found in Fig. [S] 

Tab. [J displays the canonical and quasiparticle spectra 
for the drip-line nucleus ^"Cr. In the canonical basis, [cq | 
is associated with a 3si/2 state and is about 180 keV. The 
next low- lying state (2c?5/2) is within an energy interval 
of AE « 500 keV. Those two states are separated from a 
core of orbitals by £" « 3.5 MeV. Equivalently, the sepa- 
ration energy in the quasiparticle basis is |eo | ~ 430 keV, 
whereas four quasiparticle states (si/2 and are with 
an energy spread of AE w 470 keV, and are further sep- 
arated from higher-excited states by E' « 3.2 MeV. The 
separation energy Sn for ^''Cr is compatible with the phe- 
nomenological binding energy necessary for the appear- 
ance of light halo nuclei, namely 2 MeV/A^ « 137 keV. 
According to the discussion of Sec. IIV B 21 the energy 
scales at play in the three last bound Cr isotopes corre- 
spond to ideal halo candidates. 

The criteria introduced in Sec. IIV B 41 are now applied. 
Fig. [19] shows the average number of nucleons participat- 
ing in the potential halo. Whereas A^haio is consistent 
with zero for < 50, a sudden increase is seen beyond 
the iV = 50 shell closure. The existence of a decorrelated 
region in the density of the last three Cr isotopes is con- 
sistent with the evolution of the neutron densities along 
the isotopic chain in Fig. [501 For N > 50, such a behav- 
ior translates into a non-zero value of A'haio- The value 
of -/Vhaio remains small in comparison to the total neu- 
tron number, as the decorrelated region is populated by 
~ 0.45 nucleons on the average in ^°Cr. In absolute value 
however, A^haio is comparable to what is found in light 
s-wave halo nuclei like ^^Be, where roughly 0.3 nuc leons 
constitute the decorrelated part of the density [ll5j | . 

The halo factor (Si?haio is shown in Fig. [21] as a func- 
tion of A/N. The halo contributes significantly to the 



Can. spectrum *°Cr 



Exc. spectrum '^^Cr 



El 



AEi 



E' 



3si 
2d. 



lS9/2 
1/5/2 
1/5/2 
2pi/2 



[MeV] 



> 

-0.178 
-0.670 



-4.062 
-8.676 
-8.676 
-8.942 

< -10 



/5/2 
Pl/2 
99/2 



E' 



AE < 



El 



E:: ~ A" [MeV] 



> 10 
8.694 
8.960 
4.103 



^5/2 


0.893 


^5/2 


0.832 


Sl/2 


0.728 


Sl/2 


0.427 







Table L (Color Online) Neutron canonical energies e" in *°Cr 
and separation energies = E" — A", as predicted by the 
{SLy4-|-REG-M} functional. Quasiparticle states with a spec- 
troscopic factor smaller than lO"'^ are not included. 
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Figure 19: (Color Online) Average number of nucleons partic- 
ipating in the halo along the Cr isotopic chain, as a function 
of the nuclear mass, as predicted by the {SLy4-f REG-M} 
functional. Theoretical uncertainties are included (see text). 



total neutron r.m.s. radius (up to 0.13 fm) beyond 
the V = 50 shell closure. The latter result can be recast 
as a splitting of the total r.m.s. radius into a core and a 
halo contributions, as displayed in Fig. [22] In contrast to 
the Helm model, shell effects are here properly separated 
from halo ones, e.g. the core r.m.s. radius includes a 
kink at V = 50 which is due to the filling of least bound 
states and not to the halo per se. Only the physics re- 
lated to the existence of truly decorrelated neutrons is 
extracted by A'haio and (5i?haio- The kink of the neutron 
r.m.s radi us (i) was not assumed as a halo signature a 
priori [t^, IIIGI ] but recovered a posteriori (ii) must be 
corroborated using finer tools such as A'haio and (5i?haio 
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r [fm] 

Figure 20: (Color Online) Neutron densities for even-even Cr 
isotopes, from ^''Cr to ^^Cr. The proton density of ^^Cr is 
given (dashed-dotted line) as a reference for the neutron skin. 
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Figure 22: (Color Online) Total neutron root-mean-square 
radius (solid line) and core contribution (dashed line) for 
chromium isotopes, as predicted by the {SLy4-|-REG-M} 
functional. 
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Figure 21: (Color Online) Halo factor parameter SRhaio in the 
Cr isotopic chain. 



to extract quantitatively the contribution of the halo to 
that kink. 

To characterize further this halo region, individual con- 
tributions -/Vhaio,i are evaluated. The results are summa- 
rized in Tab. |TT1 As expected, the main contributions 
to the halo come from the most weakly-bound states, 
while for non-halo nuclei, like ^''Cr, all contributions are 
consistent with zero. At the neutron drip-line, impor- 
tant contributions are found from both 3si/2 and 2d^/2 
states. The latter £ = 2 states contribute for almost 50% 
of the total number of nucleons in the decorrelated re- 
gion, although this state is more localized than the 3s 1/2 
because of its binding energy and of the effect of the cen- 
trifugal barrier. Such hindrance effects are compensated 
by the larger canonical occupation of the /2 states and 
the larger intrinsic degeneracy of the shell. The signif- 
icant contribution of the 1 = 2 states could not be ex- 
pected from the standard qualitative analysis presented 



in Sec. IIIBI or. with a few exceptions [1(1, from the ex- 
perience acquired in light nuclei. 

Finally, the probability Pi for nucleons occupying the 
canonical state 0" to be in the outer region r > tq in ^°Cr 
is typical of s— wave halo systems; i.e. 49% for the 3si/2 
state and a little bit lower for the 2^5/2 state, around 
26%. 

The analysis method applied to neutron-rich Cr iso- 
topes demonstrates unambiguously that a halo is pre- 
dicted for the last three bound isotopes. We have in- 
deed been able to characterize the existence of a decor- 
related region in the density profile for isotopes beyond 
the iV = 50 shell closure. Such a region contains a small 
fraction of neutrons which impact significantly the ex- 
tension of the nucleus. It is generated by an admixture 
of S112 and c?5/2 states, whose probabilities to be in the 
halo region r > tq are compatible with what is seen in 
light halo nuclei. This provides the picture of a rather 
collective halo building up at the neutron drip-line for Cr 
isotopes. 



C. Sn isotopes 

So far, the prediction of halos in tin isotopes beyond 
the iV = 82 shell closure [59|] have been based on the 
Helm model, whose limitations have been pointed out in 
Sec. IIII B 21 The robust analysis tools introduced in the 
present work are expected to give more reliable results. 
Neutron densities of Sn isotopes do exhibit a qualitative 
change for TV > 82, as seen in Fig. [221 However, the tran- 
sition is smoother than in the case of chromium isotopes 
(Fig. I^D)) . This is partly due the increase of collectivity 
associated with the higher mass. There are also spe- 
cific nuclear-structure features that explain the absence 
of halo in drip-line Sn isotopes. 

Tab. [HT] displays the canonical and quasiparticle spec- 
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Table IL Contributions of the least bound canonical orbitals 
to the number of nucleons in the decorrelated region, and 
probabilities for those states to be in the outer region r > ro. 
The data are provided for the four last (predicted) bound Cr 
isotopes. 



tra for the drip-linc nucleus ^^^Sn. The energy scales at 
play are not compliant with the definition of a halo, as 
can also be seen from Fig. [Ml In the canonical basis, 
the separation energy E is roughly 1.2 MeV, whereas six 
states with an energy spread AE w 3.8 MeV are sepa- 
rated from a core of orbitals by a gap E' « 5.5 MeV. 
Equivalently in the quasiparticle basis one has (i) 
Sn = E ^ 1.5 MeV, (ii) four low-lying quasiparticles with 
a spread AE « 3.4 MeV (iii) separated from higher ex- 
citations by E' ^ 5.6 MeV. The energy spread of the 
low-lying states AE is too large to favor the forma- 




Figure 23: (Color Online) Same as Fig.[20]for the Sn isotopes. 
The "separation" between the two groups of neutron densi- 
ties occurs for N = 82. Proton density of ^"^"Sn is given as a 
reference in dashed-dotted line. 



tion of a halo. Also, according to the phenomenologi- 
cal criterion extracted for light halo nuclei, the separa- 
tion energy of ^^^Sn should have been of the order of 
2 MeV/A2/3 « 64 keV for a halo to emerge. 
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Figure 24: (Color Online) Same as Fig. [18] for neutron sepa- 
ration energies of Sn isotopes. 

The iVhaio parameter is displayed in Fig.[55| The max- 
imum value of A'halo, around 0.18, is very small compared 
to the total number of nucleons. The absolute numbers 
are also smaller than the ones obtained in (lighter) Cr ha- 
los. We may add that the value of A'halo found here is of 
the same order of magnitude as those encountered for a 
non-halo p-wave nucleus such as ^^N, where around 0.12 
neutron out o f six reside in average in the classically for- 
bidden region jll5l ] . An interesting feature is the decrease 
of A^haio for A^ > 166. This is a consequence of the filling 
of the highly degenerate lii3/2 state right at the drip- 
line (see Fig. [5]). As the number of neutrons occupying 
the lii3/2 shell increases, the depth of the one-body po- 
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Table III: Same as Tab. Ufor the neutron canonical energies 
of ^^"'Sn, and associated separation energies |e"| of ^'^^Sn. 



tential also increases and the shells become more bound, 
thus more localized. As this happens over a significant 
number of neutrons, the effect on A^haio is visible. This 
constitutes an additional hindrance to the formation of 
halos from low-lying high angular-momentum states. 

The second halo parameter SRhaio displayed in Fig. 1261 
shows that the decorrelated region has little influence on 
the nuclear extension, of the order of 0.02 fm. Its con- 
tribution is found to be much less than predicted by the 
Helm model. The heavy mass of tin isotopes hinders the 
possibility of a sharp separation of core and tail contri- 
butions in the total density and thus, of the formation of 
a halo. 

The analysis of single-particle contributions, summa- 
rized in Tab. IIVI confirms the latter analysis. First, 
3pi/2i 3p3/2 and 2/7/2 {£ = 3) states contribute roughly 
the same to A'haio- For higher angular-momentum or- 
bitals, the effect of the centrifugal barrier is seen: the 
lhg/2 and lji3/2 orbitals, the latter being the least bound 
orbital, do not contribute significantly to the decorrelated 
region. Finally, individual probabilities Pi remain very 
small, and do not exceed a few percent. 




Figure 25: (Color Online) Average number of nucleons in the 
spatially decorrelated region for Sn isotopes. For comparison, 
A^haiol^^Cr) is shown as a horizontal dashed-dotted line. 
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Figure 26: (Color Online) Halo factor parameter 57?haio in the 
Sn isotopic chain. For comparison purposes, the maximum 
value of SRhaio obtained for Cr isotopes is represented as a 
horizontal dashed-dotted line. 
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Figure 27: (Color Online) Same as Fig. [22]for Sn isotopes. 
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Table IV: Same as Tab. HH for Sn isotopes. 



For all the reasons exposed above, only a neutron skin 
effect is seen in tin isotopes, and no significant halo for- 
mation is envisioned. Of course, all results presented here 
have been obtained with a particular EDF and it is of in- 
terest to probe the sensitivity of the predictions to the 
different ingredients of the method [58| . 

In any case, the two previous examples already pro- 
vide a coherent picture regarding the properties of halo or 
non-halo medium-mass nuclei. In particular, it is rather 
obvious that the notion of giant halo [tI, [tI, [tI, [tI, 



Isol . [8l| constituted of six to eight neutrons is mislead- 
ing. Indeed, such a picture was obtained by summing 
up the total occupations of loosely bound orbitals. Al- 
though loosely bound orbitals are indeed responsible for 
the formation of the halo, nucleons occupying them still 
reside mostly inside the nuclear volume. It is thus unap- 
propriate to simply sum up their occupations to charac- 
terize the halo. The identification of the halo region in 
the presently proposed method led us to define the more 
meaningful quantity A'haio- 



VI. CONCLUSIONS 

The formation of halo in finite many-fcrmion systems 
is a quantum phenomenon caused by the possibility for 
non-classical systems to expand in the classically forbid- 
den region. One difficulty to further understand this phe- 
nomenon resides in the absence of tools to characterize 
halo properties in a quantitative way. Light nuclei con- 
stitute an exception considering that the quantification 
of halo properties in terms of the dominance of a clus- 
ter configuration and of the probability of the weakly- 
bound clusters to extend beyond the classical turning 
point is well acknowledged [SJ, [H, [H, [s^j • Several at- 
tempts to characterize halos in systems constituted of 
tens of fermions have been made but were based on loose 
definitions and quantitative criteria. Such a situation 
is unsatisfactory because important questions, such as 
the very existence of halos at the neutron drip-line of 
medium-mass nuclei, are still open. 

After demonstrating the inability of the Helm model 
to provide reliable predictions, a new quantitative anal- 
ysis method has been developed to identify and charac- 
terize halos in finite many-fermion systems in a model- 
independent fashion. It is based on the decomposition of 
the internal one-body density in terms of overlap func- 
tions. The definition of the halo, as a region where nucle- 
ons are spatially decorrelated from the others, has been 
shown to be connected to specific patterns of the inter- 
nal one-body density and of the energy spectrum of the 
{N — l)-body system. In particular, halos can be char- 
acterized by the existence of a small nucleon separation 
energy E, a small energy spread AE of low-lying excita- 
tions, and a large excitation energy E' of the upper- lying 
states with respect to low-lying bunched ones, in the ex- 
citation spectrum of the {N — l)-body system. 

Based on the new analysis method, it is possible to 
extract the radius rp beyond which the halo, if it exists, 
dominates over the core. Such an identification of vq has 
been validated by extensive simulations. It is important 
to stress that the method does not rely on an a priori 
separation of the density into core and halo components. 
The latter are extracted from the analysis itself, using 
the total matter density as the only input. Several quan- 
titative observables are then introduced, namely (i) the 
average number of fermions participating in the halo, (ii) 
the influence of the halo region on the total extension 



21 



of the system, and (iii) the contributions of individual 
overlap functions to the halo. 

The new analysis method has been applied to the 
results obtained from energy density functional calcu- 
lations of chromium and tin isotopes using the code 
HFBRAD [s^l- Drip- line Cr isotopes appear as ideal halo 
candidates whereas tin isotopes do not. 

For drip- line Cr isotopes, the average fraction of nu- 
cleon participating in the halo is of the order of ~ 0.5. 
Such a v alue is compliant with those found for light halo 
systems |ll5j |. The halo region was also found to influ- 
ence significantly the nuclear extension. Contributions 
from several individual components, including i — 1 ones, 
were identified, contradicting the standard picture aris- 
ing from few-body models. The notion of collective halos 
in medium-mass nuclei has been introduced. 

In the case of Sn isotopes, the average number of nu- 
cleons participating in the halo is very small and has no 
influence on the nuclear extension. Thus, the drip-line 
phenomenon discussed previously for tin isotopes [59| is 
rather a pronounced neutron skin effect. Such skin effects 
are of course of interest as they emphasize the isovector 
dependence of the energy density functionals. However, 
they should not be confused with halo systems which dis- 
play an additional long tail of low density matter. 

This preliminary study on two isotopic series gives 
promising results and validates the theoretical grounds of 
the analysis. With upcoming new radioactive beam facil- 
ities, interaction cross-sections are expected to be mea- 
surable in the drip- line region of Z « 26 elements |ll7j . 
This would constitute a giant leap towards an extensive 
comparison between theoretical and experimental works 
on drip-line physics. 
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Appendix A: INTERNAL ONE-BODY DENSITY 

1. Definition 

In the laboratory frame, the one-body density is the 
expectation value of the operator 



which leads for the TV-body ground state to 

p(f) = iVy"dri...dfAr_iK(fi...fW_i,f)|2 
= dS^l...d^N^2dRN-l 

x|*^(6...6v-2,r-i?Ar_i)p (A2) 

where $f (fi . . . rW) = $f . . . fw-i). Using that l>f 
is invariant under translation of the system, one eas- 
ily proves that the one-body density in the laboratory 
frame is also translationally invariant, p(r-t-a) = p(r), 
and thus is uniform. This is a general property of transla- 
tionally invariant systems which underlines that the den- 
sity in the laboratory frame is not the proper tool to 
study self-bound systems. 

The relevant object for self-bound systems is the inter- 
nal one-body density matrix, defined as the expectation 
value of the operator 

N 

P[i](r,r") = 6{Rn)J2 Sir- h + Rn-i) 

i=l 

X H Sih-r'i). (A3) 

k,l = l..N 

where 

1 ^ 

(A4) 

The internal density defined with respect to the center- 
of-mass of the remaining {N — l)-body^^^ '' is of direct 
relevance to knockout reactions [104, llOSl . IllSj . Using 
the orthogonality relationship [lOoj j 

J dfi... dfN $f ... tat) S{Rn) (ri . . . fjv) = % , 

(A5) 

and dH]), one obtains [lO^ . [lO^ . [l07t 

Pli]{r,f") ^ N j dri . . .drN-i^Q* [ri . . .TN^i.r') 

x6{Rn-i)^o in ■ ■ - rN-ur) 
= Y^^UnMr), (A6) 



^ One could define another internal one-body density taking the 

p(r) = ^ ^ 5{r — Ti) , center-of-mass of the Af-body system as a pivot point. This is a 



more relevant choice to analyze electron scattering data. 
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which shows that the internal one-body density ma- 
trix is completely determined by internal overlap func- 
tions [l06l |. 

The internal one-body density /5[i] (r ) is the local part 
of the internal density matrix, and is the expectation 
value of the operator 

N 

mir) = S{RN)Y.S{f- r- + R},_,) . (A7) 
•i=i 

According to Eq. (|A6[) . one has 
Pmir) = E \^'^(^^)\" - E ^^I^^WP ■ (A8) 



2. Nuclear EDF calculations 

The behavior of the internal one-body density high- 
lighted in Sec. IIVI is general and model-independent. It 
is valid for any finite many-fermion system, as long as 
the inter-fermion interaction is negligible beyond a cer- 
tain relative distance. Of course, when an approximate 
treatment of the iV-body system is used, a certain deteri- 
oration of the properties of the density can be observed. 
In the case of EDF calculations however, some more pro- 
found issues are raised. 

First, an important clarification regarding the physi- 
cal interpretation of the quantities at play in the calcu- 
lations must be carried out. In single-reference imple- 
mentations of the nuclear EDF method, one manipulates 
the so-called "intrinsic" one body density, in the sense 
that it is built from an auxiliary state that breaks sym- 
metries of the nuclear Hamiltonian, e.g. translational, 
rotational and gauge invariance. The intrinsic density 
is associated with a wave packet from which true eigen- 
states, and their laboratory and internal densities, can be 
recovered by restoring broken symmetries through multi- 
reference EDF calculations [4l[. In practice, the intrinsic 
density is used as a good approximation to the internal 
density, e.g. when analyzing electron scattering data. 
Still, the intrinsic density of a symmetry breaking state 
and the internal density associat ed w ith the true eigen- 
state of inter est a re different (ll9j . As a result, EDF 
methods^ [ll3| |. expressed directly in term s of the in- 
ternal density are currently being considered |ll3j . 

As just said, the EDF intrinsic density has been shown 
in many cases to be a good approximation of the in- 
ternal density extracted through electron scattering. In 



In shell model, the internal wave-function is explicitly computed 
when the center-of-mass part of the A'^ body wave function can 
be mapped onto a Os state. 

The SR-EDF method, as it is currently applied to self-bound 
nuclei, is not related to an existence theorem a la Hohenberg- 
Kohn. 



practice, one identifies the lower component of the in- 
trinsic HFB wave- function VJ{f) with the internal over- 
lap function (pHr) leading from the ground state of the 
TV-body system to the corresponding excited state of the 
{N — l)-body system In particular, and this is key 
to the present discussion, the asymptotic part of VJ*(r) 
satisfies the free Schrodinger equation [82| . just as the 
asymptotic part of does. The smallest energy 

IcqI thus extracted relates to the exact separation en- 
ergy, i.e. an analogue to Koopmans' theorem derived 
originally in the case of Hartree-Fock approximation ap- 
plies. Given that the intrinsic density fEq. lM]) expressed 
in terms of the lower component of HFB quasiparticle 
wave-functions reads the same as the internal density ex- 
pressed in terms of overlap functions (Eq. lASp , the analy- 
sis method developed in Sec. IIVI including the occurrence 
of crossing patterns, applies directly to the former. 



a. Slater determinant as an auxiliary state 

In the implementation of the EDF method based on 
a Slater determinant, explicit spectroscopic factors are 
either zero or one, and behave according to a step func- 
tion = 6(ep — el). The single-particle orbitals (pi are 
identified with overlap functions and the density takes 
the form given by Eq. (|A8p . 



b. Quasiparticle vacuum as an auxiliary state 

In the implementation of the EDF method based on a 
quasiparticle vacuum, the one-body density can be eval- 
uated using either the canonical states 4>1 or the lower 
components V^' of the quasiparticle states 

i V 

(A9) 

where j'' relates to the total angular momentum. In the 
present case, the spectroscopic factor identifies with 
the quasiparticle occupation Nj^ defined by Eq. ([7]). This 
underlines that implementation of the EDF approach 



It can bo shown that the porturbative one-quasiparticle state 
•r]i^^) contains N + Ui^ — fi^ particles on the average if |<I>) is 
constrained to A'^ particles on the average. It is only for deep- 
hole quasiparticle excitations {vf « 1) that the final state will be 
a good approximation of the (A^ — l)-body system. The correct 
procedure, that also contains some of the rearrangement terms 
alluded to above, consists of constructing each one-quasiparticle 
state self-consistently by breaking time-reversal invariance and 
requiring (A^ — 1) particles in average, or of creating the quasi- 
particle excitation on top of a fully paired vacuum designed 
such that the final state has the right average particle num- 
ber [1201 . Il2ll | . The overlap functions and spectroscopic factors 
can be computed explicitly in such a context. 
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based on a quasiparticle vacuum incorporates exp licitl y 
parts of the spreading of the single-particle strength jl22l |. 
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Figure 28: (Color Online) Neutron quasiparticle occupation 
N" as a function of the separation energy in *°Cr, calcu- 
lated with the {SLy4-|-REG-M} functional. Conventions from 
Fig. [2] are used to label individual quasiparticle states. Only 
quasiparticles with occupations greater than 10"'^ are dis- 
played. 



The function S*" — /(|e"|), whose typical behavior is 
presented in Fig. [55] for ^°Cr, takes values between zero 
and one. The difference between hole-like quasiparticle 
excitations and particle-like ones is visible. Indeed, S*" 
increases with excitation energy |e"| for hole- like exci- 
tations. This constitutes the main branch which tends 
towards a step function when correlations are not explic- 
itly included into the auxiliary state; i.e. for the EDF ap- 
proach based on an auxiliary Slater determinant. Spec- 
troscopic factors of particle-like quasiparticle excitations 
remain small and go to zero for high-lying excitations. 



Appendix B: DETERMINATION OF THE HALO 
REGION 



Let us start with a very crude toy model, where every- 
thing is analytical. The total density p is assumed to be 
a superposition of a core pc and a tail p^, both taking 
the form 



(Bl) 



This amounts to considering that the asymptotic 
regime is reached in the region of the crossing between 
Pc and ph, and we neglect for now the factor. In this 
model the second-order (base-10) log-derivative of the to- 
tal density is analytical, as well as the exact positions of 
(i) its maximum rmax (h) the point rg where the halo 
density is exactly equal to ten times the core one. Then, 
the ratio 7^(^o) = log"g /9(ro)/ log'/g p(rmax) can be eval- 
uated and becomes in the weak binding limit of interest 

Kh/Kc 



(B2) 



This shows that the position where there is a factor of 
ten between pc and ph is equivalently obtained by finding 
the position where there is a given ratio between the value 
of the second-order log-derivative of the density and its 
maximal value. The critical value 40/121 « 0.33 found 
in the toy model is not believed to be accurate for com- 
plex nuclei, as (i) the asymptotic regime is not reached 
at the crossing point and is more complicated because 
of the r^^ factor (ii) the total density is a superposition 
of more than two components. However, we expect the 
one-to-one correspondence between ratios on the densi- 
ties and ratios on log'/p p to hold in realistic cases. Thus, 
the position where the halo dominates the core by one or- 
der of magnitude can be found using log'/p p as the only 
input. 

More realistic model calculations have been used to 
characterize the position of ro- The total density is taken 
as a linear combination of core and halo contributions. 
Their relative normalization are free parameters in this 
simulation, allowing to artificially change the fraction of 
halo in the total density 



Ptot{r) = NcPc{r) 



(B3) 



where and Nh ~ -/V^ are the number of nucleons in 

the core part and in the halo part, respectively. The den- 
sities Pc and pi, are normalized to one. We considered (i) 
simple models, where the core and each halo components 
are defined as 



1 



r < Ro, 



(B4) 



r> R 



: 



Mi standing for a normalization constant. This model 
only accounts for the basic features of the nuclear den- 
sity: a uniform core of radius Rq and a spatial exten- 
sion becoming larger as a; — s- (ii) double Fermi models, 
where the un-physical sharp edge in the logarithmic rep- 
resentation of the previous density is smoothened out 



Po 



(B5) 



1 



(iii) semi-phenomenological models, which fulfill the 
asymptotic behavior of Eq. ((28)) . Core and tail densities 
vanish at r = 0, as well as their derivatives with respect 
to X, y and z, in order to avoid singularities at the nu- 
cleus center 



z, m order to avoid smguiarities at the nu- 
123. In Refs. [l23 . [T25lJT26l] . such densities 
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were adjusted on experimental data. The core part was 
defined as 



Pc{r) = 



1 + 





a. 


2 





, (B6) 



where a = 1 for neutrons, and the halo density as 

^2 



Ph{r) = Po,t 



e "t , 



(B7) 



(iv) more realistic models, where the core density is still 
defined as in Eq. (jB6p . but the halo contributions are 
realistic wave functions taken from self-consistent EDF 
calculations of Cr and Sn isotopes. 
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Figure 29: (Color Online) [Main panel] Ratio between the 
second-order log-density at ro and its peak value log'/o p(rmax) 
[Top panels] Distribution of r/ro for which TZir) is equal 
to a given value (left panel: 7^(r) = 0.35, middle panel: 
7^(r) = 0.4, right panel: 7^(r) = 0.5), and ro. 

The results from a wide panel of test cases are pre- 
sented in Fig. [221 We recall that ro is defined as the 



radius for which Phiro) = 10pc{ro). The goal is to de- 
termine such a radius through the value of another ratio: 
TZ{ro). For each simulation, the position tq and the ratio 
7?.(ro) are computed. The main panel of Fig. shows the 
distribution of Tl{ro) which is peaked around 0.4. This 
value is greater than in the toy model case for the reasons 
detailed above. The distribution of 7?.(ro) is asymmetric, 
but the tail towards high values corresponds precisely to 
non-halo systems. In any case, we are going to reflect 
such an asymmetry into the choice of theoretical uncer- 
tainties in the determination of rg. 

The inserts of Fig. [29] display the distribution of ratio 
r/ro corresponding to a given value of 7?,(r). The ra- 
tio TZ{r) = 2/5 (top-center panel) indeed picks out quite 
consistently the radius ro. For TZ{r) — 40/121 (top-left 
panel), the position r is in most cases below ro. As a con- 
sequence, the average ratio between tail and core compo- 
nents in the density will be consistently below ten in this 
case. On the contrary for TZ{r) = 1/2 (top-right panel), r 
is systematically larger than ro, meaning that the Ph/Pc 
larger than ten on the average. In the end, it appears 
that ro is indeed well picked out through the condition 



7^(r•o) 



(B8) 



For those reasons, we use those values of 72. (r) to set 
the error bars on the determination of ro. 

Of course, we need to account for the fact that a dif- 
ference by one order of magnitude between core and halo 
densities to define the halo region is somewhat arbitrary 
and that the corresponding radius ro cannot be perfectly 
picked out in all cases through Eq. (jBSp . As a result, we 
add a tolerance margin to the definition of ro by allow- 
ing 7^(ro) to vary between 40/121 « 0.35 and 1/2. The 
upper margin is greater than the lower one to account 
for the asymmetry of the peak in Fig. [2S1 Note that the 
procedure chosen to determine ro combined with that 
asymmetry put us on the safe side, i.e. the radius found 
through that procedure, if not perfect, is likely to be too 
large, leading to a slight underestimation of the halo fac- 
tors Nb_aio and i5i?haio- 
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